Synonymous sites for accessibility around microRNA binding sites in bacterial spot and speck disease resistance genes of tomato.

The major causes of mass tomato infections in both covered and open ground are agents of bacterial spot and bacterial speck diseases. MicroRNAs (miRNAs) are 16-21 nucleotides in length, non-coding RNAs that inhibit translation and trigger mRNA degradation. MiRNAs play a significant part in plant resistance to abiotic and biotic stresses by mediating gene regulation via post-transcriptional RNA silencing. In this study, we analyzed a collection of bacterial resistance genes of tomato and their binding sites for tomato miRNAs and Pseudomonas syringe pv. tomato miRNAs. Our study found that two genes, bacterial spot disease resistance gene (Bs4) and bacterial speck disease resistance gene (Prf), have a 7mer-m8 perfect seed match with miRNAs. Bs4 was targeted by one tomato miRNA (sly-miR9470-3p) and three Pseudomonas syringe pv. tomato miRNAs (PSTJ4_3p_27246, PSTJ4_3p_27246, and PSTJ4_3p_27246). Again, Prf gene was found to be targeted by two tomato miRNAs namely, sly-miR9469-5p and sly-miR9474-3p. The accessibility of the miRNA-target site and its flanking regions and the relationship between relative synonymous codon usage and tRNAs were compared. Strong access to miRNA targeting regions and decreased rate of translations suggested that miRNAs might be efficient in binding to their particular targets. We also found the existence of rare codons, which suggests that it could enhance miRNA targeting even more. The codon usage pattern analysis of the two genes revealed that both were AT-rich (Bs4 = 63.2%; Prf = 60.8%). We found a low codon usage bias in both genes, suggesting that selective restriction might regulate them. The silencing property of miRNAs would allow researchers to discover the involvement of plant miRNAs in pathogen invasion. However, the efficient validation of direct targets of miRNAs is an urgent need that might be highly beneficial in enhancing plant resistance to multiple pathogenic diseases.

Turning off a few overexpressed genes in prostate cancer with microRNAs using a 7mer-seed match model.

PURPOSE: This research aims to identify the miRNAs that could target the genes overexpressed in prostate cancer so that miRNA-based therapeutics could be developed. METHODS: A 7mer-m8 model of microRNA targeting was utilized in order to analyse the relationship between microRNAs and overexpressed genes. The efficiency of miRNA binding was investigated using various parameters namely free energy (AMFE), GC and GC3 content, translation efficiency, cosine similarity metric, mRNA stability, free energy of RNA duplex, and base compositional difference. BLAST2GO software was used to elucidate the functional roles of the genes overexpressed in prostate cancer. RESULTS: The current research reveals that the coding sequences of the genes were found targeted with multiple miRNAs. For instance, the HPN gene was targeted by the microRNA miR-4279 at two distinct sites i.e. 263-278 and 746-761 in the coding sequence. In the present study, it was observed that the target region of the genes exhibited a comparatively high GC and GC3 contents in comparison to the flanking regions. A low translational rate and weak relationship between RSCU and tRNA were obtained which may be due to the absence of optimal codons. CONCLUSION: In this study, we have uncovered the human miRNAs that have potential for binding to the coding sequences of 14 most overexpressed genes in prostate cancer and thereby could silence those genes.

Comparison of compositional constraints: Nuclear genome vs plasmid genome of Pseudomonas syringae pv. tomato DC3000.

Pseudomonas syringae pathovar (pv) tomato DC3000 (PsPto) causes bacterial speck disease in Solanum lycopersicum (tomato) and Arabidopsis thaliana. This bacterium has attracted a lot of attention from researchers. Codon usage bias (CUB) analysis is a helpful tool to understand the gene and genome biology of an organism and also provides insights into the evolutionary mechanism. Several studies have been conducted on different bacteria, but seldom on PsPto. In the present study, the codon usage pattern of three transcriptomes of PsPto, i.e., nuclear genome, plasmid A and plasmid B were systematically analyzed. The PsPto transcriptomes were found to be strongly biased towards G and C nucleotides at the third codon position. ENC, neutrality plots, PR2 plot and COA analysis also suggested that the CUB of PsPto is mainly affected by natural selection. A comparative analysis based on the relative codon deoptimization index (RCDI) for adaptation of plasmid A and plasmid B to the nuclear genome revealed that plasmid B was closer to the nuclear genome.

Methane sources from waste and natural gas sectors detected in Pune, India, by concentration and isotopic analysis.

Methane (CH4) is a potent greenhouse gas and also plays a significant role in tropospheric chemistry. High-frequency (sub-hourly) measurements of CH4 and carbon isotopic ratio (δ13CH4) have been conducted at Pune (18.53°N, 73.80°E), an urban environment in India, during 2018-2020. High CH4 concentrations were observed, with a mean of 2100 ± 196 ppb (1844-2749 ppb), relative to marine background concentrations. The δ13CH4 varied between -45.11 and -50.03 ‰ for the entire study period with an average of -47.41 ± 0.94 ‰. The diurnal variability of CH4 typically showed maximum values in the morning (08:00-09:00 local time) and minimum in the afternoon (15:00 local time). The deepest diurnal amplitude of ~500 ppb was observed during winter (December-February), which was reduced to less than half, ~200 ppb, during the summer (March-May). CH4 concentration at Pune showed a strong seasonality (470 ppb), much higher than that at Mauna Loa, Hawaii. On the other hand, δ13CH4 records did not show distinct seasonality at Pune. The δ13CH4 values revealed that the significant sources of CH4 in Pune were from the waste sector (enhanced during the monsoon season; signature of depleted δ13CH4), followed by the natural gas sector with a signature of enriched δ13CH4. Our analysis of Covid-19 lockdown (April to May 2020) effect on the CH4 variability showed no signal in the CH4 variability; however, the isotopic analysis indicated a transient shift in the CH4 source to the waste sector (early summer of 2020).

Identification of promising CD8 and CD4 T cell epitopes for peptide vaccine formulation against SARS-CoV-2.

The novel virus "severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)" has been responsible for the worldwide pandemic causing huge devastation and deaths since December 2019. The disease caused by this virus is known as COVID-19. The present study is based on immuno-informatics approach to develop a multi-epitope-loaded peptide vaccine to combat the COVID-19 menace. Here, we have reported the 9-mer CD8 T cell epitopes and 15-mer CD4 T cell epitopes, free from glycosylation sites, belonging to three proteins, viz. surface glycoprotein, membrane glycoprotein and envelope protein of this virus. Immunogenicity, aliphatic amino acid, antigenicity and hydrophilicity scores of the predicted epitopes were estimated. In addition, other physicochemical parameters, namely net charge, Boman index and amino acid contents, were also accounted. Out of all the epitopes, three CD8 T cell epitopes viz. PDPSKPSKR, DPSKPSKRS and QTQTNSPRR and three CD4 T cell epitopes viz. ASYQTQTNSPRRARS, RIGNYKLNTDHSSSS and RYRIGNYKLNTDHSS were found to be efficient targets for raising immunity in human against this virus. With the help of our identified potent epitopes, various pharma industries might initiate efforts to incorporate those epitopes with carrier protein or adjuvant to develop a multi-epitope-loaded peptide vaccine against SARS-CoV-2. The peptide vaccines are usually cost-effective and therefore, could be administered as a preventive measure to combat the spread of this disease. Proper clinical trials must be conducted prior to the use of identified epitopes as vaccine candidates.

A Study on microRNAs Targeting the Genes Overexpressed in Lung Cancer and their Codon Usage Patterns.

As reported by WHO in 2018, there were 2.09 million victims of lung cancer and 1.76 million fatalities worldwide. Tobacco remains the biggest hazard in causing this lethal disease. To execute the computational analysis, the overexpressed lung cancer genes were retrieved from literature and subsequently their complete coding sequences (CDS) were downloaded. The mature microRNA sequences of human were extracted from miRBASE. The 7mer-m8 perfect seed match between miRNAs and mRNAs was found. Following filtration, 7 genes were selected that possessed binding sites for maximum miRNAs viz., MUC5B (miR-4479, miR-1227-5p, miR-3940-5p, miR-604, miR-4455, miR-4267, miR-6750-3p, miR-4530, miR-5587-5p, miR-4508, miR-4534, miR-4443, miR-4253, miR-1321, miR-4655-5p, miR-4297, miR-4296, miR-1268a, miR-3178, miR-4750-3p, miR-1306-3p, miR-1268b, miR-3656, miR-1233-3p, miR-6804-5p), MUC16 (miR-4456, miR-1205, miR-665, miR-6808-3p, miR-1279, miR-4257, miR-1227-5p, miR-888-3p, miR-4455, miR-4267, miR-4294, miR-1275, miR-4288, miR-1178-5p, miR-4314, miR-6829-3p, miR-548av-5p, miR-1294, miR-5587-5p, miR-3622b-5p, miR-1273f, miR-4770, miR-4327, miR-4318, miR-4531, miR-4534, miR-4443, miR-7106-5p, miR-3125, miR-3650, miR-4325, miR-4266, miR-7976, miR-1290, miR-4500, miR-7160-5p, miR-4291, miR-1306-3p, miR-6130, miR-4430, miR-4725-5p, miR-4441, miR-6077, miR-1304-5p, miR-7515, miR-3182, miR-6134), COL1A1 (miR-3665, miR-1227-5p, miR-6132, miR-2861, miR-4530, miR-3155b, miR-3155a, miR-1292-3p, miR-4497), COL5A1 (miR-7162-5p, miR-3665, miR-6809-3p, miR-4313, miR-4531, miR-4532, miR-3155b, miR-4323, miR-1207-3p, miR-4260, miR-6071, miR-4710, miR-7162-5p), CELSR2 (miR-7150, miR-4308, miR-6132, miR-4770, miR-4534, miR-4492, miR-3960, miR-3178, miR-4291, miR-563), COL7A1 (miR-665, miR-6730-3p, miR-1227-5p, miR-4265, miR-6829-3p, miR-4297, miR-4532, miR-3181, miR-4310, miR-4441, miR-4497, miR-1237-3p), and FAT2 (miR-4267, miR-1275, miR-4770, miR-1825, miR-6895-5p, miR-4535, miR-4493, miR-940, miR-6861-3p, miR-4310, miR-4710, miR-4447, miR-4472). The miRNA-target site and their flank regions were compared with respect to site accessibility, translational rate, and relationship between RSCU and tRNAs. Higher accessibilities to miRNA-binding regions and lower translational rates indicated that miRNAs' binding to their respective targets might be efficient. The presence of rare codons might further augment miRNA targeting. The codon usage bias study of the genes related to lung cancer revealed non-uniform usage of nucleotides and comparatively higher GC content. Lower biasness prevailed in the genes and selective constraint mostly governed them. Lastly, the functionalities of target genes were also revealed. The silencing characteristic of miRNAs might be exploited to design miRNA-mediated therapy that might potentially repress the overexpressed genes in carcinoma.

Binding sites of miRNA on the overexpressed genes of oral cancer using 7mer-seed match.

The microRNAs having a length of ~ 19-22 nucleotides are the small, non-coding RNAs. The evolution of microRNAs in many disorders may hold the key to tackle complex challenges. Oral cancer belongs to the group of head and neck cancer. It occurs in the mouth region that appears as an ulcer. In this study, we collected information on the overexpressed genes of oral cancer. The coding sequences of the genes were derived from NCBI and the entire set of human microRNAs present in miRBASE 21 was retrieved. The human microRNAs that can target the overexpressed genes of oral cancer were determined with the aid of our in-house software. The interaction between microRNAs and the overexpressed genes was evaluated with 7mer-m8 model of microRNA targeting. The genes DKK1 and APLN paired with only one miRNA i.e., miR-447 and miR-6087, respectively. But the genes INHBA and MMP1 were found to be targeted by 2 miRNAs, while the genes FN1, FAP, TGFPI, COL4A1, COL4A2, and LOXL2 were found to be targeted by 16, 5, 9, 18, 29, and 11 miRNAs. Subsequently, several measures such as free energy, translation efficiency, and cosine similarity metric were used to estimate the binding process. It was found that the target region's stability was higher than the upstream and downstream zones. The overexpressed genes' GC contents were calculated, revealing that the codons in target miRNA region were overall GC rich as well as GC3 rich. Lastly, gene ontology was performed to better understand each gene's involvement in biological processes, molecular function, and cellular component. Our study showed the role of microRNAs in gene repression, which could possibly aid in the prognosis and diagnosis of oral cancer.

Influencing elements of codon usage bias in Birnaviridae and its evolutionary analysis.

Birnaviridae is a family of double stranded (ds) RNA virus with non-enveloped virions and 2-segmented genomes. These viruses are known to cause diseases in many hosts. Virus of this family has affected the fish and poultry economy in a wide sector. Unevenness in the use of synonymous codons for a particular amino acid in the coding strand of DNA is known as codon usage bias (CUB). Codons that code the same amino acid are used with variable frequency in a variety of life forms. To understand the pattern of CUB in Birnaviridae, we carried out bioinformatics study to understand the properties of coding sequences of proteins. ENC value of Birnaviridae suggested low CUB. Nucleotide analysis revealed high GC content. Parameters such as RSCU values, nucleotide skewness, translational selection, parity plot and neutrality plot were studied to investigate the pattern of codon use and it was clear that both mutational pressure and natural selection contributed to the designing of CUB in Birnaviridae family. The neutrality plot revealed natural selection to dominate the structuring of CUB and hence remained the major CUB determinant in Birnaviridae. Outcome of our study exemplified the pattern of codon use in the Birnaviridae genomes and contributed the basic primary data for fundamental evolutionary research on them.

Human tooth enamel carbon and oxygen stable isotope dataset from chalcolithic Inamgaon (India).

The tooth enamel from the human remains of ten archaeological individuals belonging to a chalcolithic site at Inamgaon, District Pune, Maharashtra, were analysed for stable carbon and oxygen isotope compositions. The human remains of the involved individuals come from three consecutive periods: Period I (1600-1400 BC; n = 2), Period II (1400-1000 BC; n = 4), and Period III (1000-700 BC; n = 4). Enamel carbonate of twenty teeth (n = 20), two from each individual, were analysed to understand the inter- and intra-individual variations in isotope ratios across the three habitational periods. The acquired dataset will help in understanding isotope baseline values for the region in the prehistoric context. The subsequent research works in the region can reuse our data in collation with other datasets for comparative investigations.

Identification of putative microRNAs in the complete genome of Mycobacterium avium and their possible interaction with human transcripts.

The grievous adversity regarding Mycobacterium avium is its ubiquitous nature. Isolation of the bacteria from drinking water, house dust, and soil, etc., is an alarming issue for the scientific community. The microRNAs are the molecular influencers of gene expression that act during the process of post transcription. A few reports claimed the existence of microRNAs or microRNA-like molecules in the prokaryotic species. Biogenesis of bacterial miRNAs requires their transport into the host cell. Subsequently, the host-encoded enzymes are exerted for the formation of bacterial mature miRNAs and their regulation. In our study, the screening of complete genome of Mycobacterium avium revealed six putative precursor microRNA sequences bearing typical secondary structures. The mature microRNAs were predicted in both arms of the secondary structures. A total of 12 possible mature microRNAs were identified in this study. The likely targets of the predicted mature miRNAs were searched in human 3' UTR. In the human transcriptome, 193 genes were possibly targeted by 12 mature miRNAs of Mycobacterium avium. The essential functionalities of the target genes included signal transduction, immune system, DNA binding, and response to stress.

Compositional features and pattern of codon usage for mitochondrial CO genes among reptiles.

The phenomenon of non-random occurrence of synonymous nucleotide triplets (codons) in the coding sequences of genes is the codon usage bias (CUB). In this study, we used bioinformatic tool kit to analyze the compositional pattern and CUB of mitogenes namely COI, COII and COIII across different orders of reptiles. Estimation of overall base composition in the protein-coding sequences of COI, COII and COIII genes of the reptilian orders revealed an uneven usage of nucleotides. The overall count of A nucleotide was found to be the highest while the overall count of G nucleotide was the least. The CO genes across the three reptilian orders were prominently AT biased. Comparison of the GC proportion at each codon position displayed that GC1 percentage ranked the highest in all the three CO genes of the reptilian orders. SCUO values indicated weaker CUB, while considerable variation of SCUO values existed in the three CO genes across the studied reptiles. Relative synonymous codon usage (RSCU) values indicated that mostly the A ending codons were preferred. Based on the parameters namely neutrality plot, mutational responsive index and translational selection, we could conclude that natural selection was the major evolutionary force in COI, COII and COIII genes in the studied reptilian orders. However, correspondence analysis, parity plot and correlation studies indicated the existence of mutation pressure as well on the CO genes.

Identification of novel microRNAs in Rous sarcoma Virus (RSV) and their target sites in tumor suppressor genes of chicken.

A small non-coding, evolutionarily conserved regulatory RNA molecule known as microRNA (miRNA) regulates various cellular activities and pathways. MicroRNAs remain evolutionarily conserved in different species of same taxa. They are present in all organisms including viruses. Viral miRNAs are small, less conserved and less stable and have higher negative minimal folding free energy than miRNAs of different organisms. The size of viral precursor miRNA is approximately 60-119 nucleotides in length. The structure of the mature miRNA sequences is predicted by using higher negative MFE (ΔG) value. Rous sarcoma Virus (RSV), named after its inventor Peyton Rous, has been known for causing tumors in the chicken for which it is known as an oncogenic retrovirus. Using specific criteria we have predicted 5 potential miRNAs in RSV which targeted 8 tumor suppressor genes in Gallus gallus. This study aims to predict the potential miRNAs, secondary structures and their targets for better understanding of the regulatory network of Rous sarcoma virus miRNA in forming sarcoma.

Extreme rainfall deficits were not the cause of recurring colonial era famines of southern Indian semi-arid regions.

Using information contained in the eighteenth to twentieth century British administrative documents, preserved in the National Archives of India (NAI), we present a 218-year (1729-1947 AD) record of socioeconomic disruptions and human impacts (famines) associated with 'rain failures' that affected the semi-arid regions (SARs) of southern India. By mapping the southern Indian famine record onto long-term spatiotemporal measures of regional rainfall variability, we demonstrate that the SARs of southern India repeatedly experienced famines when annual rainfall reduced by ~ one standard deviation (1 SD), or more, from long-term averages. In other words, 'rain failures' listed in the colonial documents as causes of extreme socioeconomic disruptions, food shortages and human distress (famines) in the southern Indian SARs were fluctuations in precipitation well within the normal range of regional rainfall variability and not extreme rainfall deficits (≥ 3 SD). Our study demonstrates that extreme climate events were not necessary conditions for extreme socioeconomic disruptions and human impacts rendered by the colonial era famines in peninsular India. Based on our findings, we suggest that climate change risk assessement should consider the potential impacts of more frequent low-level anomalies (e.g. 1 SD) in drought prone semi-arid regions.

Analysis of codon usage bias in mitochondrial CO gene among platyhelminthes.

The phenomenon of non-uniform usage of the synonymous codons, where some codons are given more preference to others, is known as codon usage bias (CUB). CUB is known to be determined by two major evolutionary forces i.e. mutation pressure and selection. We used various approaches to understand the codon usage pattern in mitochondrial CO (MT-CO) genes involved in complex IV of the respiratory chain (RC) as no work was reported yet. Our present study revealed that CUB was relatively high and the coding sequences were rich in A and T. Correspondence analysis further indicated that A/T compositional properties under mutational pressure might be affecting the codon usage pattern and was different in different classes for MT-CO gene. A highly significant correlation between A% and A3%, T% and T3%, G% and G3%, C% and C3%, GC% and GC3% in all the classes indicated that compositional constraints under mutational pressure and natural selection might affect the CUB. Neutrality plot indicated that both natural selection and mutational bias affected the CUB, where, natural selection played the major role as compared to mutational pressure.

Analysis of codon usage of Horseshoe Bat Hepatitis B virus and its host.

In the present analysis, codon usage strategies and base distribution of Horseshoe bat hepatitis B virus (HBHBV) were analyzed and compared with its host Rhinolophus sinicus, as no work was yet reported. The magnitude of synonymous codon usage bias (CUB) in the virus and its host was low with higher proportion of the base C. Notably, 21 more frequently used codons, 19 less frequently used codons and 3 underrepresented codons (TCG, ACG and GCG) were found to be similar in both virus and its host coding sequences. Neutrality plot analysis reported greater role of natural selection in HBHBV (67.84%) and R. sinicus (76.90%) over mutation pressure. Base skewness and protein properties also influenced the CUB of genes. Further, codon usage analysis depicted, HBHBV and R. sinicus had many similarities in codon usage patterns that might reflect viral adaptation to its host.

Free energy of mRNA positively correlates with GC content in chloroplast transcriptomes of edible legumes.

In the present study, the results of nucleotide composition analysis showed that the legume chloroplast (cp) transcriptomes were AT rich. From the neutrality plot, we observed that natural selection might have played a major role, while mutation pressure played a minor role in the CUB of cp transcriptomes. Highly significant (p < 0.05) negative correlation was found between mRNA free energy (mFE) and scaled chi-square for entire mRNA in Cicer arietinum and Lens culinaris suggesting that the release of higher energy by entire mRNA molecule might be associated with higher degree of codon usage bias in these two crop plants. Further, highly significant (p < 0.01, p < 0.05) positive correlation of mFE for entire mRNA was found with GC3 and that of mFE for 39 bases with GC, GC1, GC2 and GC3 contents among all the legumes. This indicated that higher GC content might induce the release of more free energy by cp transcriptomes.

Silencing lung cancer genes using miRNAs identified by 7mer-seed matching.

Lung cancer (LC) is the main cause of cancer-associated deaths in both men and women globally with a very high mortality rate. The microRNAs (miRNAs) are a class of noncoding RNAs consisting of 18-25 nucleotides. They inhibit translation of protein through binding to complementary target mRNAs. The non-coding miRNAs are recognized as potent biomarkers for detection, development and treatment of malignancy. In this study, we screened a set of 12 genes over expressed in small cell lung cancer, non small cell lung cancer and the genes involved in both categories and their binding sites for human miRNAs as no work was reported yet. Screening of human miRNAs revealed that a few genes showed numerous miRNA binding sites. Free energy values of mRNA sequences revealed that they might acquire compact folded structure causing complexity for miRNAs to interact. GC content in the target site was relatively higher than that of their flanks. It was observed through analysis of cosine similarity metric and compAI parameters that the genes related to lung cancer were encoded with non optimal codons and thus might be translationally less efficient for producing polypeptides. Gene ontology analysis was carried out to understand the diverse functions of these 12 genes.

Diurnal and seasonal variability of CO2 and CH4 concentration in a semi-urban environment of western India.

Amongst all the anthropogenically produced greenhouse gases (GHGs), carbon dioxide (CO2) and methane (CH4) are the most important, owing to their maximum contribution to the net radiative forcing of the Earth. India is undergoing rapid economic development, where fossil fuel emissions have increased drastically in the last three decades. Apart from the anthropogenic activities, the GHGs dynamics in India are governed by the biospheric process and monsoon circulation; however, these aspects are not well addressed yet. Towards this, we have measured CO2 and CH4 concentration at Sinhagad, located on the Western Ghats in peninsular India. The average concentrations of CO2 and CH4 observed during the study period are 406.05 ± 6.36 and 1.97 ± 0.07 ppm (µ ± 1σ), respectively. They also exhibit significant seasonal variabilities at this site. CH4 (CO2) attains its minimum concentration during monsoon (post-monsoon), whereas CO2 (CH4) reaches its maximum concentration during pre-monsoon (post-monsoon). CO2 poses significant diurnal variations in monsoon and post-monsoon. However, CH4 exhibits a dual-peak like pattern in pre-monsoon. The study suggests that the GHG dynamics in the western region of India are significantly influenced by monsoon circulation, especially during the summer season.

Composition, codon usage pattern, protein properties, and influencing factors in the genomes of members of the family Anelloviridae.

The present study was carried out on 62 genome sequences of members of the family Anelloviridae, as there have been no reports of genome analysis of these DNA viruses using a bioinformatics approach. The genes were found to be rich in AC content with low codon usage bias (CUB). Relative synonymous codon usage (RSCU) values identified the preferred codons for each amino acid in the family. The codon AGA was overrepresented, while the codons TCG, TTG, CGG, CGT, ACG, GCG and GAT were underrepresented in all of the genomes. A significant correlation was found between the effective number of codons (ENC) and base constraints, indicating that compositional properties might have influenced the CUB. A highly significant correlation was observed between the overall base content and the base content at the third codon position, indicating that mutations might have affected the CUB. A highly significant positive correlation was observed between GC12 and GC3 (r = 0.904, p < 0.01), which indicated that directional mutation pressure influenced all three codon positions. A neutrality plot revealed that the contribution of mutation and natural selection in determining the CUB was 58.6% and 41.4%, respectively.

Understanding the codon usage patterns of mitochondrial CO genes among Amphibians.

A universal phenomenon of using synonymous codons unequally in coding sequences known as codon usage bias (CUB) is observed in all forms of life. Mutation and natural selection drive CUB in many species but the relative role of evolutionary forces varies across species, genes and genomes. We studied the CUB in mitochondrial (mt) CO genes from three orders of Amphibia using bioinformatics approach as no work was reported yet. We observed that CUB of mt CO genes of Amphibians was weak across different orders. Order Caudata had higher CUB followed by Gymnophiona and Anura for all genes and CUB also varied across genes. Nucleotide composition analysis showed that CO genes were AT-rich. The AT content in Caudata was higher than that in Gymnophiona while Anura showed the least content. Multiple investigations namely nucleotide composition, correspondence analysis, parity plot analysis showed that the interplay of mutation pressure and natural selection caused CUB in these genes. Neutrality plot suggested the involvement of natural selection was more than the mutation pressure. The contribution of natural selection was higher in Anura than Gymnophiona and the lowest in Caudata. The codons CGA, TGA, AAA were found to be highly favoured by nature across all genes and orders.